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Abstract 

In analogy with parton distribution functions, also parton fragmen- 
tation functions obey matching conditions when crossing heavy-flavour 
thresholds. We compute these matching conditions at next-to- leading 
order in the strong coupling constant a s in the MS scheme. Our re- 
sults can be used for the dynamical generation of the heavy-flavour 
component in next-to-leading order fits to light-hadrons fragmentation 
functions. Furthermore, when computing perturbatively the charm 
fragmentation function from first principles and evolving it to higher 
scales, our matching conditions should be used for consistency when 
crossing the bottom threshold. 



1 Introduction 



In full analogy with structure functions, also fragmentation functions change 
when heavy-flavour thresholds are crossed. In a system with n L light flavours, 
and a heavy flavour of mass m, the inclusive production of a hadron at a scale 
much below the heavy-flavour threshold is given by 

£- = [dzdp T J2 Dt\z^)^t^5(p T - zp T ) , (1) 
dp T J ^ dp T 

where I nL is the set of all light partons 

In L = {qi,qi, ■■■q nL ,Qn Il , 9} , ( 2 ) 

dal nh) /dp T are the short distance cross sections, calculable in perturbation 
theory, and Df 1 ^ are the partonic fragmentation functions that describe the 
hadronization of the light partons % into the hadron. The partonic fragmen- 
tation functions D^ 1 ^ are normally not calculable perturbatively, but their 
evolution in terms of the scale ft is given by the Altarelli-Parisi equations 
with n h flavours. 

Renormalization is performed in the so called decoupling scheme [1]: di- 
vergences associated with light particles are subtracted in the MS scheme, 
while divergent heavy-quark loop effects are subtracted at zero momentum. 
This results in the fact that, for momenta much below the heavy-quark mass, 
the heavy-quark existence can be completely ignored. In the following, we 
will also refer to this scheme as the n L scheme. 

Much above the heavy-flavour threshold, an equation of the form (0) also 
holds, but with n h replaced by n — n L + 1. The obey Altarelli-Parisi 
equations for n flavours, and the standard MS scheme is used for all flavours. 
We will refer to this scheme as the n scheme. 

For a heavy enough quark, the relation of D^ n ' with Z)( nL ) can be com- 
puted in perturbative QCD. At the leading logarithmic level, we have 

D (n) ( s _ f Dt\z,Li) + 0(a s ) for i G I„ L , . 

1 /i) " \ OK) for i G {h, h} Wh6n ~ m ' (3) 

where we have introduced the notation h = q n for the heavy flavour. The 
matching condition Q simply states that the presence of the heavy quark 
has effects of order a s in processes at scales near its mass. A relation of the 
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same form as (j2J) holds for parton densities below and above a heavy-flavour 
threshold. In this case, however, the next-to-leading order (NLO) accurate 
matching condition has also been known for a long time [2]. It requires 
to be of order a% and the difference of and D^ nL ^ for the light flavours 
to be at most of order for \i exactly equal to m. This matching condition 
is commonly used in next-to-leading logarithmic (NLL) parton distribution 
functions fits, where the heavy-quark distributions c(x, ji) and b(x, /x) are 
radiatively generated, rather than directly fitted to the data. 

For fragmentation functions, NLO matching has never appeared in the 
literature, and it has never been used in global fits. Usually [3-8] the heavy- 
flavour fragmentation functions are simply parametrized and fitted to the 
data. In [9], they are instead generated dynamically, but using only the 
leading-order matching condition ©. 

In the present work, we compute the matching conditions for fragmen- 
tation functions at the NLO level in the MS scheme. The NLO matching 
conditions should be used for consistency in NLL fragmentation-function fits 
that generate heavy flavours dynamically. In particular, in the computation 
of the charm fragmentation function from first principles [10], the match- 
ing conditions should be used when crossing the bottom threshold. In this 
framework, at low energies (i.e. not much above the charm mass), the charm 
is treated as a heavy quark, in order to provide a perturbative expression for 
its fragmentation function. Near the bottom threshold, the bottom is treated 
as heavy, while all other quarks (including charm) are considered light. 



2 Matching conditions at NLO accuracy 

We consider a light-hadron H inclusive production process in e + e~ annihila- 
tion 1 

e+e" -> H + X (4) 

at a center-of-mass energy Q much above the heavy-flavour mass, so that 
power-suppressed m/Q terms can be neglected, but not too large, so that 
a s \og{m/Q) < 1. 

In this regime, we can perform the calculation in both the decoupling 
and the full MS scheme. In the first scheme, we can neglect to resum large 

1 Universality of the fragmentation function guarantees that the result obtained in a 
specific process remains valid for all processes. 
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logarithms of m/Q, since as^og(m/Q) is small, use the fragmentation 
functions, and include the heavy-flavour effects by fixed 0(a s ) calculation. 
On the other hand, in the second scheme, we can neglect power-suppressed 
m/Q terms, treat all flavours as light, and use the full MS scheme, with the 
£)( n ) fragmentation functions and the corresponding factorization formulae. 

In the decoupling scheme, we can write the differential cross section for 
the production of a light hadron H with energy fraction x = 2Eh/Q as 



(5) 

The first term involves the fragmentation function and the short distance 
cross sections <7j for all partons i, excluding the heavy flavour. The second 
term is the 0(a s ) heavy-flavour contribution, arising from direct production 
of a heavy-flavoured pair, followed by gluon radiation and fragmentation, as 
depicted in Fig. The cross section can be easily obtained from the 




Figure V. Heavy-flavour contribution to the production of a light hadron. 
appendix of Ref. [11] 

— T — = °hh 7T- C F 2 <^ log — + log(l - y) - 1 \ , (6) 

ay Z7T 2/1 m J 

where a h h is the Born cross section for the production of a heavy quark- 
antiquark pair, y = 2E g /Q is the energy fraction carried by the gluon, and 
C F = 4/3. Terms suppressed by powers of m/Q have been consistently 
neglected in Eq. ©. 

The calculation for the same observable da/dx can also be performed us- 
ing the full (massless) MS scheme, treating all flavours as massless. Defining 
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l n = I nL U {h, h}, we obtain 



da 

dx 



Uy 



EA (n) (^,M) 



dy 



iein 



(7) 



The difference between the expressions (JHJ) and (J2J) yields 



'dy 



+ 



V 

'dy 

y 



Wix/y^^-D^ix/y^) 



Df\x/y^) ^^ -D^\x/y^) 



dy 

d&i nh) (y,») 



+ - 



x dy 
y 



i£{h,h} 



dy 



D^\x/y,fi) 



dy 

dvhhgjy) 
dy 



0.(8) 



In the first term of Eq. (JHJ) we have replaced <7> and af 1 ^ with <7j, since 
these cross sections only differ by terms of order a%. In the second term, 
since a g is of order a s , we can consider Dg™^ = D g n \ since their difference, 
according to Eq. (J3J), is of order a s . Furthermore 



da^\y,fi) dfr^\y,n) _ da hhg (y,ii) 



dy 



dy 



dy 



(9) 



where da h f ig (y, n)/dy is the massless MS-subtracted cross section for the pro- 
cess of Fig. ^ i- e - a gluon emitted from the (now only nominally) heavy 
quark. 

Finally, in the third term of Eq. (JHJ), are of order a s (according to 
Eq. (JHJ)), so that we only need the Born term for the hard cross section 



dy 



a hh 6(l-y)+0(a s ) . 



(10) 



Since Eq. (JHJ) must hold for arbitrary electric charges of each quark flavour, 
it follows immediately that, neglecting terms of order a|, 



Df l \x, fj) = D™{x, fj) for i e I nL , i ^ g 



(11) 
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We thus find 



—D g (x/y,fi) 



x y 



da h h g {y,n) da hh Jy) 



dy 



+ / i E ^W) 

J x y 



which implies 



D { :\x,f,) + D^ (x,fM) 



(«)/ 



iE{h,h} 



dy 

a hh Jx V 



hhg \ 

dy 
dy 



(12) 



d(r h - hq (y) da hh Jy,n) 



' hhg 

dy 



dy 



The MS subtracted cross section d'hhgi.U ^ 1^) is given by [IX, 12] 



D g (x/y,/i) . 

(13) 



1 = ^ TT" i 2 logy + log(l - y) + log — 



and thus, using Eqs. (jOJ and (JHJ), we obtain 

1 + (1-Z/) 2 



(14) 



Z^W) = 4 n W) 



/ —D g (x/y,n)^C F 

Jx y 2tt 



?7 



log— ? - 1 - 2 logy 



.(15) 



Equations ()11)1 and (jl5j) are the matching conditions for the D q and .D^, 
accurate up to terms of order a s . We remark that the right hand side in 
Eq. (J15|) does not vanish for \i — m. Under this respect, matching condi- 
tions for parton fragmentation functions differ from those for distribution 
functions [2], which instead vanish at the heavy-quark threshold. 

Contrary to D q and D^, the matching condition for the gluon fragmenta- 
tion function D g has not yet been determined to NLO accuracy. In fact, in 
e + e~ annihilation, we are only sensitive to the leading order value of D g , since 
gluon production is suppressed by a power of a s . On the other hand, if we 
consider a wide evolution span from m to Q, such that a s log(Q/m) « 1, an 
0(a s ) error in D g at m would propagate into an 0(a s ) error in all the other 
components of D through evolution. We must therefore provide a matching 



condition for D g accurate to order a s . We will show that, in fact, 



D^{z^)=D ( Q n -\z^) + 0{al) for/, 



m 



(16) 



5 




which is the same result that holds for parton densities [2]. In order to 
prove Eq. we consider the process of light-hadron production from a 

gluon. We can imagine that the gluon is produced in some physical process. 
For concreteness, we can think of a "super-heavy" quark produced in e + e~ 
annihilation (renormalized in the same way in the n L and n scheme), that 
radiates a gluon. The gluon produces a light hadron by fragmentation, ac- 
cording to the graph a in Fig. |21 We now want to consider the difference 

b 
d 



Figure 2: Contributions to the production of a light hadron from a gluon. 

in the radiative corrections to this process in the n L and in the n scheme. 
It is clear that processes where the gluon splits into a pair of light partons 
give the same contribution in both schemes. These are in fact processes of 
order a 3 , diagrammatically identical in the n L and n scheme. In these pro- 
cesses (that carry an extra power of a s ) the strong coupling constant and the 
fragmentation functions can be taken at leading order, and therefore they 
coincide, at 0(a s ), in the two schemes. The process g — > hh followed by the 
fragmentation of the heavy quark (shown in graph b of the figure) exists only 
in the n scheme, but it is suppressed by a^, since the fragmentation function 
Dh is of order a s , and therefore does not give a relevant contribution. The 
only corrections that can contribute to the difference are shown in graphs c 
and d. In c, a heavy quark loop corrects the gluon propagator. We include 
in c also the corresponding renormalization counterterm. Graph d represents 
the difference in the collinear subtractions in the n and n L schemes. 

Calling II the contribution from the heavy-quark loop in graph c, and C r 
the renormalization counterterm, we have 

n (n L ) + C (n L ) = q ; ( 17 ) 
H (n) + C (n) = C (n) = y ^ 

Equation (J 17)) follows from the fact that the heavy-quark loop is subtracted 
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at zero momentum, and the gluon is on shell in c. Equation (|18|) follows from 
the fact that a massless quark loop at zero momentum is zero in dimensional 
regularization, so only the standard counterterm for a quark loop in the MS 
scheme survives. 

For graph d, the difference in the collinear subtraction counterterm arises 
only from the P gg term of the Altarelli-Parisi splitting functions, which is 
the only term that depends explicitly upon the number of light flavours. We 
have 

-is^'M-WlH^r- (19) 

that cancels exactly the remnant from Eq. (jl8j) . The cancellation is in fact 
obvious. The fermion contribution of the vacuum polarization from a quark 
loop vanishes, and thus its ultra-violet (UV) component must cancel exactly 
its infra-red (IR) component. The renormalization counterterm should be 
equal to the UV component with the opposite sign, and thus equals the 
IR component. Finally, the collinear counterterm arises precisely from the 
IR component of a quark loop, and so it compensates the renormalization 
counterterm. 

The only remaining contribution that could differ in the n and n L schemes 
is the one represented in graph a. It is proportional to the product of the 
gluon fragmentation function times the coupling a s arising from the emission 
of the gluon. Since the total result must be scheme independent, we must 
therefore have 

4 n V)^ n) (^) = 4 nL V)^ nL) (^) • (20) 

From the matching condition for the running coupling a s 

a\> = a K s 1 + — — log — , 21 
\ 37r m 2 / 

we therefore infer 

Df\z^) = D^\ Z ^){l-^lo g ^ . (22) 

Summarizing, for fi w m and up to corrections of order a\ we have 

D h)r^ x ^) = [ — D 9( x /y^) 
j x y 
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a s l + (l-y) 2 
2vr y 



log— ^ - 1 - 2 logy 



i/i 



i/i 



for i = q u ...,q. 



For future reference, we provide the Mellin transform of -D^l. Defining 



h/K 



(23) 

(24) 
(25) 

(26) 



we have 

l(n) 



iv(^-i) l log ^" 1) + 

2 



(JV-1): 



AT2 



+ 



(iV + 1) 2 



D 9 (N) 



(27) 



3 Conclusions 

In this paper we have calculated, at order a s and in the MS scheme, the 
matching conditions for parton fragmentation functions at the heavy-quark 
thresholds. Such conditions can be used to generate radiatively and to next- 
to-leading accuracy charm and bottom contributions to light-hadron or pho- 
ton production via fragmentation. This is achieved by evolving light-quarks 
and gluon fragmentation functions from a low- to a high-energy scale, and 
through the heavy-quark thresholds. 

In a similar fashion, these matching conditions at next-to-leading ac- 
curacy are required for full consistency when evolving the perturbatively- 
calculated charm fragmentation function through the bottom threshold. 

The possibility of generating dynamically the heavy-quark contributions 
to fragmentation will allow to fit light-hadron fragmentation functions with 
less free parameters than presently done, in full analogy with modern parton 
distribution functions analyses. We defer these fits to a future paper. 
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